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Abstract 

Following our investigation of the USA Standard and Poor index anti-bubble that 
started in August 2000 [Quantitative Finance 2, 468-481 (2002)], we analyze thirty 
eight world stock market indices and identify 21 "bearish anti-bubbles" and 6 
"bullish anti-bubbles." An "anti-bubble" is defined as a self-reinforcing price trajec- 
tory with self-similar expanding log-periodic oscillations. Mathematically, a "bearish 
antibubble" is characterize by a power law decrease of the price (or of the logarithm 
of the price) as a function of time and by expanding log-periodic oscillations. We 
propose that bearish anti-bubbles are created by positive price-to-price feedbacks 
feeding overall pessimism and negative market sentiment further strengthened by 
inter-personal interactions. Bullish anti-bubbles are here identified for the first time. 
The most striking discovery is that the majority of European and Western stock 
market indices as well as other stock indices exhibit practically the same log-periodic 
power law anti-bubble structure as found for the USA S&PSOO index. These anti- 
bubbles are found to start approximately at the same time, August 2000, in all 
these markets. This shows a remarkable degree of synchronization worldwide. The 
descent of the worldwide stock markets since 2000 is thus an international event, 
suggesting the strengthening of globalization. 
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1 Introduction 



Financial bubbles are loosely defined as phases of over-valuations of stock 
market prices above the fundamental prices. Such over-valuations may be in 
accord with the theory of rational expectations, leading to the concept of 
rational expectation bubbles [2,3,23,25,42], may be due to exogenous causes 
or sunspots (see for instance [5]) or may result from a variety of departures 
from pure and perfect agent rationality [32,33,31,34,36]. 

A series of papers based on analogies with statistical physics models have pro- 
posed that most financial crashes are the climax of the so-called log-periodic 
power law signatures (LPPS) associated with speculative bubbles resulting 
from imitation between investors and their herding behavior [39,16,22,13,40]. 
In addition, a large body of empirical evidence supporting this proposition 
have been presented [41,39,13,19,17,40]. A complementary line of research has 
established that, while the vast majority of drawdowns occurring on the ma- 
jor financial markets have a distribution which is well- described by an ex- 
ponential or a slightly fatter distribution in the class of stretched exponen- 
tials, the largest drawdowns are occurring with a significantly larger rate than 
predicted by extrapolating the bulk of the distribution and should thus be 
considered as outliers [14,40,20,12]. A recent work [21] has merged these two 
lines of research in a systematic way to offer a classification of crashes as 
either events of an endogenous origin associated with preceding speculative 
bubbles or as events of an exogenous origin associated with the markets re- 
sponse to external shocks. Two hallmarks of criticality have been documented: 
(i) super-exponential power law acceleration of the price towards a "critical" 
time tc corresponding to the end of the speculative bubble and (ii) log-periodic 
modulations accelerating according to a geometric series signaling a discrete 
hierarchy of time scales. Globally over all the markets analyzed, Ref. [21] iden- 
tified 49 outliers, of which 25 were classified as endogenous, 22 as exogenous 
and 2 as associated with the Japanese "anti-bubble" . Restricting to the world 
market indices, Ref. [21] found 31 outliers, of which 19 are endogenous, 10 
are exogenous and 2 are associated with the Japanese anti-bubble. The exoge- 
nous crashes, not preceded by LPPS could be in each case associated with an 
important piece of information impacting the market. 

All these results taken together formulate a general hypothesis according to 
which imitation between investors and their herding behavior lead to specu- 
lative bubbles of financial markets with accelerating overvaluation decorated 
by accelerating oscillatory structures possibly followed by crashes or change of 
regimes. The key concept is the existence of positive price-to-price feedbacks. 
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When speculative prices go up, creating wealth for some investors, this may 
attract other investors by word-of-mouth interactions, fuelling further price 
increases. This in turn promotes a wide-spread interest in the media which 
promotes and amphfies the self-fulfiUing wishful thinking [29] , with seemingly 
reasonable or rational theories advanced to justify the price increases. These 
processes generate more investor demand, fuelling further the expansion of 
the speculative bubble. The positive price-to-price feedback mechanism has 
recently been formulated mathematically in a nonlinear generalization of the 
Black-Scholes stochastic differential equation [37] and in a nonlinear model of 
stock market prices combining the positive price-to-price feedback with nonlin- 
ear negative feedback due to fundamental trading together with inertia [11,38]. 
It was there shown that the speculative bubble becomes unstable, reflecting 
the fact that high prices are ultimately not sustainable, since they are high 
only because of expectations of further price increases. The bubble eventually 
bursts, and prices come faUing down. The feedback that fed the bubble carries 
the seeds of its own destruction, and so the end of the bubble and the crash 
are often unrelated to any really significant news on fundamentals [36]. 

The same feedback mechanism may also produce a "negative" bubble or "bear- 
ish anti-bubble," that is, downward price movements propelling further down- 
ward price movements, enhancing pessimism by inter-personal interactions. 
Johansen and Sornette [15] proposed indeed that such imitation and herding 
mechanism may also lead to so-called "anti-bubbles" with decelerating mar- 
ket devaluations following market peaks. The concept of "anti-bubble" was 
introduced to describe the long-term depression of the Japanese index, the 
Nikkei, that has decreased along a downward path marked by a succession of 
ups and downs since its all-time high of 30 Dec. 1989 [15,18]. The concept 
of anti-bubble restores a certain degree of symmetry between the speculative 
behavior of the "bull" and "bear" market regimes. This degree of symmetry, 
after the critical time tc, corresponds to the existence of "anti-bubbles," char- 
acterized by a power law decrease of the price (or of the logarithm of the 
price) as a function of time t > tc, down from a maximum at tc (which is 
the beginning of the anti-bubble) and by decelerating/expanding log-periodic 
oscillations [15,18]. Another anti-bubble was found to describe the gold fu- 
ture prices after its all-time high in 1980. The Russian market prior to and 
after its speculative peak in 1997 also constitutes a remarkable example where 
both bubble and anti-bubble structures appear simultaneously for the same 
tc- Several other examples have been described in emergent markets [17]. 

In a recent paper [43], we have uncovered a remarkable similarity in the be- 
havior of the US S&P500 index from 1996 to August 2002 and of the Japanese 
Nikkei index from 1985 to 1992 (11 years shift), with particular emphasis on 
the structure of the bearish phase which is qualified as an anti-bubble accord- 
ing to the previous classification. Specifically, we found the existence of a clear 
signature of herding in the decay of the S&P500 index since August 2000 with 
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high statistical significance, in the form of strong log-periodic components 
decorating a power law relaxation. 



Here, we show that the (bearish) anti-bubble that started around August 
2000 on the USA stock market is actually a world-wide phenomenon with a 
high degree of correlation and synchronization between most of the western 
markets. To our knowledge, only during the crash of October 1987 and in 
its aftermath did stock markets worldwide exhibited a similar or stronger 
correlation [1,30]. 



2 Identification of anti-bubbles in world stock market indexes 

2.1 Qualification of an anti-bubble 

Following the philosophy of Ref. [21] and references therein (see also [36] for 
a general review and references therein), we qualify an anti-bubble by the 
existence of a regime of stock market prices well-fitted by the expression 



which embodies the log-periodic power law signature. Note that the phase (f) 
does nothing but provide a time scale T since a; In (r) + = In (r/T) with 
the definition cf) — —to In T. cj) thus disappears by the choice of T as the time 
unit. This stresses the fact that, if the phase is not a fundamental parameter 
of the fit since it can be get rid of by a suitable gauge choice, it contains 
nevertheless an important information on the existence of a characteristic time 
scale, expression (1) obeys the symmetry of discrete scale invariance [35], that 
has been proposed to be a hallmark of cooperative behavior of interacting 
agents [22,13,36]. The meaning of the adjective "well-fitted" will be clarified 
below, first by presenting visual evidences in figures and then more formally 
by statistical tests. For a speculative bubble, we have 



which is the time to the end of the bubble occurring at tc- For an anti-bubble, 
we have 



which is the time since the beginning of the anti-bubbles at tc- The exponent 
a should be positive in order for the price p{t) to remain finite at t — tc 



Inp {t)^A + Bt" + Ct" cos [a; In (r) + 0] , 



(1) 



r ^tc-t , 



(2) 



T^t-t 



(3) 



4 



[13]. In general, speculative bubbles exhibiting the LPPS described by (1) 
with (2) are followed by crashes or strong corrections that are "outliers" [21]. 
It may happen that some of these speculative bubbles transform rather into 
anti- bubbles described by (1) with (3). As we said, this occurred for instance 
for the Russian speculative bubble ending in 1997. In contrast, anti-bubbles 
correspond in general to enduring corrections of stock markets that follow 
a period of strong growth, as exemplified by the trajectory of the Japanese 
Nikkei index [15,43] which culminated in Dec. 1989 and then has suffered a 
non-stop decay decorated by oscillations [15,18]. 

Following our previous finding [43] of a strong influence of a log-periodic har- 
monic at the angular log-frequency 2uo for the S&P500 index, we also present 
fits including the effect of a harmonic at 2a;. In this goal, we postulate the 
formula 

\np{t)p^A + Br"" + Cr'^ cos [uj In (r) + ^i] + Dt" cos [2uj In (r) + ^2] ,(4) 

which differs from equation (1) by the addition of the last term proportional 
to the amplitude D. The two phases 0i and 02 now define two time scales 
Ti = e^'^i/'^ and T2 = e^'^'^l'^ . It is thus no longer possible to make them 
disappear simultaneously by a choice of time unit. 

2.2 Methodology 

For each stock market described below, we use equations (1) and (4) to fit the 
logarithm of the stock market indices over an interval starting from a time 
tstart and ending in September, 30, 2002. If we knew the critical time tc, then 
an obvious choice would be tstart = tc- This choice would be optimal since it 
allows us to use expression (1) for the longer possible time span compatible 
with the occurrence of the anti-bubble. Not knowing tc precisely, in order to 
be consistent with the meaning of expression (1) with (3), we should ensure 

that tstart ^ tc- 

Furthermore, in accordance with the intuitive meaning of an anti-bubble, we 
would like to take tstart close to the last strong maximum in 2000 and then 
carry out a sensitivity analysis with respect to tstart- Fortunately, we shall show 
that the critical times tc estimated from the fit of the data for different tgtart 
do not disperse much. In other words, the fits with different tgtart are robust 
and tc is not very sensitive to igtart- We shall come back to this point later in 
Sec. 2.6 which will be focused on the best possible characterization of tc- 

As part of the sensitivity analysis with respect to tstart, we shall also use the 
following trick, which ensures that the impact of tgtart is minimized. Since 
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nothing informs us a priori about the time ordering of igtart and tc, following 
[43], we modify expression (3) into 

T=\t-t,\. (5) 



This definition (5) has the advantage of removing the constraint of tc < tstart 
in the optimization. We see this as an advantage because this constraint has 
not deep meaning and does not contain any specific information on the data 
as it results solely from the analyzing procedure and the arbitrary choices for 
tstart- It thus enables us to test the robustness of fits by scanning different igtart 
[43]. Not knowing tc precisely, in order to be consistent with the meaning of 
expression (1) with (3), we ensure that tstart > tc by trial and error: a given 
chosen tstart is accepted only if the fit gives a critical time tc < tstart- 

While the definition (3) together with the logarithmic as well as power law 
singularities associated with formula (1) imposes that tc < tstart for an anti- 
bubble, the definition (5) allows for the critical time tc to lie anywhere within 
the time series. In that case, the part of the time series for t < tg corre- 
sponds to an accelerating "bubble" phase while the part t > tc corresponds 
to a decelerating "anti-bubble" phase. Definition (5) has thus the advantage 
of introducing a degree of flexibility in the search space for tc without much 
additional cost. In particular, it allows us to avoid a thorough scanning of tgtart 
since the value of tc obtained with this procedure is automatically adjusted 
without constraint. See Ref. [43] for a discussion of the advantages and po- 
tential problems associated with this procedure using (5). Here, we use the 
definition (5) because it has proved to provide significantly better and more 
stable fits with little need to vary tstart- 

For a given tgtart; we estimate the parameters tc, a, cu, 0, A, B and C of (1) by 
minimizing the sum of the squared residues between the fit function (1) and the 
logarithm of the real index data. Following [22,13], the three linear parameters 
A, B and C are slaved to the other parameters by solving analytically a system 
of three linear equations and we are left with optimizing four free parameters. 
To obtain the global optimization solution, we employ the taboo search [6] to 
determine an "elite list" of solutions as the initial conditions of the ensuing line 
search procedures in conjunction with a quasi-Newton method. The best fit 
thus obtained is regarded to be globally optimized. A similar procedure is used 
to fit the index with formula (4). This formula has two additional parameters 
compared with (1), the amplitude D of the harmonic and its phase 02- We 
follow a fit procedure which is an adaptation of the slaving method of [22,13]. 
This allows us to slave the four parameters A, B, C and D to the other 
parameters in the search for the best fit. With this approach, we find that the 
search of the optimal parameters is very stable and provides fits of very good 
quality in spite of the remaining five free parameters. 
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We apply this procedure to 38 stock indices all over the world including: 

• eight indices in Americas (Argentina, Brazil, Chile, USA-Dow Jones, Mex- 
ico, USA-NASDAQ, Peru, and Veneruela), 

• fourteen indices in Asia/Pacific (Australia, China, Hong Kong, India, In- 
donesia, Japan, Malaysia, New Zealand, Pakistan, Philippines, South Korea, 
Sri Lanka, Taiwan and Thailand), 

• fourteen indices in Europe (Austria, Belgium, Czech, Denmark, France, Ger- 
many, Netherlands, Norway, Russia, Slovakia, Spain, Switzerland, Turkey 
and United Kingdom) and 

• two indices in Africa/Middle East (Egypt and Israel). 

By the obvious criterion to obtain at least a solution in the fitting procedure, 
we find no evidence of an anti-bubble in the following eleven indices: Austria, 
Chile, China, Egypt, Malaysia, New Zealand, Pakistan, Philippines, Slovakia, 
Sri Lanka and Venezuela. Interestingly, except for Austria and New Zealand, 
they are emergent markets in developing countries. 

The rest of the paper is thus devoted to the study of the remaining 27 indexes 
out of our initial list of 38. 

2.3 Bearish anti-bubbles 

Anti-bubbles are identified in 21 stock market indices: Netherlands, France, 
USA Dow Jones, USA NASDAQ, Japan, Belgium, Denmark, Germany, Nor- 
way, Spain, Switzerland, United Kingdom, Israel, Brazil, Hong Kong, India, 
Peru, Taiwan, Czech, Argentina and Turkey. We refer to these anti-bubble 
by the term "bearish" to stress that they are fitted by an overall decreasing 
power law (since B < and a > 0). Figures 1-21 present the fits of these 21 
bearish anti-bubbles by expression (1) and (4). The corresponding parameters 
are hsted in Tables 1 and 2. Similar bearish anti- bubbles were observed before 
in Latin- American, Asian and Western Stock markets [15,17]. We also show 
the extrapolation of these fits by formulas (1) and (4) until mid-2004, in the 
spirit of the analysis presented for the USA S&P500 index [43]. 

These figures show that the log-periodic structures arc very prominent. For 
instance, four to five log-periodic oscillations can be identified for most of the 
cases presented. However, for some indices, the log-periodic oscillations close 
to tc are strongly affected by noise for some indices and are less clear-cut. 

The majority of predicted critical times tc for the launch of the anti-bubbles 
fall between August and November, 2000: eight in August, two in September, 
three in October and one in November. This is in agreement with the deter- 
mination tc — Aug-09-2000 for the S&P500 index [43] , suggesting a worldwide 
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synchronization of the start of a bearish anti-bubble phase. 

This analysis, paralleling that presented for the USA S&P500 index [43], sug- 
gests that many of the stock markets shown here are in a phase of recovery 
that started close to the last date, September 30, 2002, used to perform the fit. 
This recovery is predicted by the extrapolations of (1) and (4) to extend until 
some time in 2003 depending upon the markets (see the figures) before a reces- 
sion resumes for a while. However, we do not claim that these extrapolations 
should be vahd beyond roughly the end of 2003. The statement apphes to the 
markets of The Netherlands, Prance, USA (Dow Jones and NASDAQ), Bel- 
gium, Denmark, Germany, Norway, Spain, Israel, Peru, and Turkey. We note 
also that there is sometimes a substantial difference between the timing of 
the recovery and the following recession predicted by expression (1) compared 
with formula (4). This is the case for the stock markets of the Netherlands, 
Prance, the USA Dow Jones, Norway, Spain, and the United Kingdom. In 
these cases, one should be careful in interpreting these extrapolations as reli- 
able forecasts. In some cases, such as for the United Kingdom and Brazil, the 
two extrapolations are so inconsistent as being meaningless. Thus, our mes- 
sage here is not so much the forecasts but instead the remarkable consistency 
in the log-periodicity of these anti-bubble phases, shown by their common 
starting dates and similar structures quantified by the power law exponents 
and the angular log-frequencies. 

In contrast, the markets of Japan, Switzerland, Hong Kong, India, Taiwan, 
Czech and Argentina are extrapolated to continue their overall descent roughly 
till the middle of 2003 or beyond before a recovery sets in. The clear log- 
periodic structure since August 2000 shown in figure 5 for the Japanese Nikkei 
index is especially interesting because this structure follows the large scale 
anti-bubble log-periodic pattern that started in January 1990 [15] and con- 
tinued at least until the beginning of 2000 [18]. A possible interpretation of 
the novel structure identified in figure 5 is that it is a sub-structure within a 
hierarchy of log-periodic patterns occurring at many different scales, as found 
for instance in Weierstrass functions (see [8] for an interpretation of Weier- 
strass functions and their generalizations in terms of log-periodicity at many 
different scales) and suggested in [7] . We can expect more generally that sim- 
ilar multiscale log-periodicity should exist in other markets. However, these 
fine structures especially at the smaller scales (three months, monthly, weekly, 
intraday, etc.) are greatly effected or even spoiled by the intrinsically noisy na- 
ture of stock market prices, due to the fact that many more effects contribute 
potentially at small scales to scramble possible signals. Only at the large time 
scales studied here can the cooperative behavior of investors be systematically 
observed. 

Another important observation is that the log-periodic oscillations and power 
law decays are distinctly different with a smaller number of oscillations and 
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much larger "noise" for Brazil, Hong-Kong, India, Peru, Taiwan, Czech repub- 
lic, Argentina and Turkey compared with the others. This may be explained 
by the presence of stronger idiosyncratic influences, such as local crises in 
South America. For these markets, the two extrapolations obtained from ex- 
pressions (1) and (4) diverge rapidly away from each other, making them quite 
unreliable. 

As shown in Table 1, the power law exponents a of the indices of Belgium 
and Argentina are significantly larger than 1, while those of Netherland, USA 
Dow Jones, Germany, Norway, Switzerland and United Kingdom are close to or 
slightly greater than 1. In absence of the log-periodic oscillations, this would 
mean that the overall shape of these indices would be concave (downward 
plunging) rather than convex (upward curvature) as they would be for < 
a <\. Large values of a > 1 implies a steep downward overall acceleration of 
the index. But in all cases when this occurs, this is compensated by a large 
amplitude of the log-periodic oscillations. In contrast, for < a < 1, the index 
initially drops fast in the early times of the anti-bubble and then decelerates 
and approaches a constant level at long times. 

To quantify the significance level of the log-periodic oscillations in these 21 
anti-bubbles, we adopt the Lomb analysis [26] on the residuals of the logarithm 
of the indices by removing the power law [17]: 

. (6) 



where r is defined in (3). If the log-periodic formula is a correct representation 
of these indices, r{t) should be a pure cosine as a function of Inr. Thus, a 
spectral analysis of r{i) as a function of the variable Inr should be a strong 
power peak. Figure 22 presents the corresponding Lomb periodograms for all 
21 indices described in table 1 and shown in figures 1-21. Most of the Lomb 
spectral peaks give a very significant signal of the existence of log-periodic 
structures [44]. 

In addition, an harmonic of a fundamental angular log-frequency uj is visible 
at 2c<j in the Lomb periodogram for many of the markets, as found previously 
for the S&P 500 index [43]. This is the justification for including an harmonic 
log-periodic oscillatory term according to (4). Table 2 lists the corresponding 
parameters of the fits of the 21 stock market indices with expression (4) and 
shows that using of formula (4) reduces the r.m.s. errors strongly for most 
of the indexes. Only for Israel is the improvement of the fit ambiguous. In 
section 2.5, we shall come back to this issue and provide rigorous and objective 
statistical tests on the relevance of log-periodicity with a single angular log- 
frequency Lo and with the addition of its harmonics at 2a;. 
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2-4 Bullish anti-bubbles 



Figures 23-28 present the fits of six stock market indices (Austrafia, Mexico, 
Indonesia, South Korea, Thailand and Russia) using formula (1) and (4), with 
the corresponding parameters hsted in Table 3. We have separated these 6 
markets from the 21 previous ones because their fits with formulas (1) and (4) 
give a positive coefficient B, corresponding to an overall increasing market at 
large time scales. We thus call them "buUish" to describe this overall increasing 
pattern. We keep the terminology "anti-bubble" to refer to the fact that the 
log-periodic oscillations are decelerating. To the best of our knowledge, the 
identification of such bullish anti-bubbles is performed here for the first time. 
Notice that there are five log-periodic oscillations for Australia, Indonesia and 
Thailand, four for and Mexico and approximately three for South Korea and 
Russia as can be seen in the figures 23-28. This means that the log-periodic 
structures in these stock market indices are quite significant and convincing 
[44]. Table 3 shows that the predicted critical times tc of the start of these 
bullish anti-bubbles are again between August and November, 2000. 

We have also performed a fit of these 6 stock market indices with expression (4) 
which accounts for the possible presence of an harmonic log-periodic oscillatory 
term. The fits are plotted in Figs. 23-28 as dashed lines, whose parameters are 
presented in Table 4. We find that the improvement of the fits using expression 
(4) compared with (1) is very significant for Australia, Korea, Indonesia and 
Thailand. 

2. 5 Statistical test of the log-periodic term 

Since expression (4) contains formula (1) as the special case D = 0, we can 
use Wilk's theorem [28] and the statistical methodology of nested hypotheses 
to assess whether the hypothesis that D — can be rejected. Similarly, we 
can also test if C = can be rejected in (1). We consider the following three 
hypotheses. 

(1) Hq-. C — 0, corresponding to use a simple and pure power law to fit the 
stock market indices; 

(2) Hi. D = 0, corresponding to the log-periodic function (1) without any 

harmonics; 

(3) H2: -D 7^ 0, corresponding to the log-periodic function (4) which includes 
an harmonics at 2co. 

Our tests presented below show that Hq can be rejected with certainty in favor 
of Hi for all the indexes which have been analyzed. Moreover, we find that Hi 
can be rejected in favor of H2 with high statistical significance for all except 
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one index (rejection level of 10^''%). We stress that Wilk's methodology of 
nested hypothesis testing automatically takes into account the competition 
between (i) the improved fit obtained by adding fitting parameters and the 
"cost in parsimony" of adding these parameters. 

The method proceeds as follows (see [40,43] for recent implementations in 
similar contexts). Assuming a Gaussian distribution of observational errors 
(residuals) at each data point, the maximum likelihood estimation of the pa- 
rameters amounts exactly to the minimization of the sum of the square over 
all data points (of number n) of the differences Sj{i) between the mathemat- 
ical formula and the data [26]. The standard deviation aj for hypothesis Hj 
with j = 0, 1, 2 of the fits to the data associated with (1) and (4) is given by 
1/n times the sum of the squares over all data points of the differences Sj°\i) 
between the mathematical formula and the data, estimated for the optimal pa- 
rameters of the fit. The log-likehhoods corresponding to the three hypotheses 
are thus given by 

Lj — —n\nV27r — nlnaj — n/2 , (7) 

where the third term results from the product of Gaussians in the likelihood, 
which is of the form 

n 

cx ncxp[-(5f)(z))V2<7|] =exp[-n/2] , 

i=l 

from the definition (j| = (1/n) Y,^=i[Sj°\i)]'^ ■ Then, according to Wilk's theo- 
rem of nested hypotheses, the log-likelihood-ratio 

r,-,+i = -2{Lj - Lj+,) = 2n{\naj - lna,+i) , j = 0, 1 (8) 

is a chi-square variable with k degrees of freedom, where k is the number of 
restricted parameters [10]. In the present case, we have k — 1. 

The Wilk test thus amounts to calculating the probability Pjj+i that the ob- 
tained value of Tjj^i can be overpassed by chance alone. If this probability 
Pj,j+i is small, this means that chance is not a convincing explanation for the 
large value of Tjj+i which becomes meaningful. This implies a rejection of the 
hypothesis that C = (resp. D = 0) is sufficient to explain the data and favor 
the fit with C 7^ (resp. D 7^ 0) as statistically significant. In other words, if 
the observed value of the probabihty 1 — Pj,j+i that Tq^i (respectively of Ti_2) 
does not exceed some high-confidence level (say, the 99% confidence level) of 
the x^, we then reject the hypothesis Hi (respectively H2) in favor of the 
hypothesis Hq (respectively Hi), considering the additional term C (respec- 
tively D) redundant. Otherwise, we accept the hypothesis Hi (respectively 
H2, considering the description with Hq (respectively Hi) insufficient. 
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For each stock index, we fit the corresponding time series starting from tgtart 
and ending on September 30, 2002 to a simple power law, to the log-periodic 
function (1) and to the formula (4) respectively, and thus obtain ctq, cti and 
(72. Then we can calculate Tjj^i from (8) and the corresponding probabilities 
Pj,j+i- The results of the Wilk tests are presented in Table 5. The values To,! 
are extremely large for all indices, which reject with extremely high statistical 
significance the hypothesis that a pure power law is sufficient compared to 
a log-periodic power law. For the test of H2 against Hi, Ti 2 is found very 
large (> 40) for most of the indices, except for Israel (Ti_2 = 2.7) and Mexico 
(Ti 2 = 25.9). Even in the case of Mexico, the improvement obtained by adding 
a the harmonic term (hypothesis H2) is nevertheless very significant since 
Ti 2 = 25.9 corresponds to a probability of rejection of H2 Pi,2 less than 
10^'*%. Thus, only for Israel, we find that D = can not be rejected at the 
confidence level of 95%. This refiect the fact that the reduction of the r.m.s. 
errors when going from formula (1) to (4) is less that 1%. The lack of significant 
improvement can also be seen visually in Fig. 13. 

Since the assumption of Gaussian noise is most probably an under-estimation 
of the real distribution of noise amplitudes, the very significant improvement 
in the quahty of the fit brought by the use of both formulas (1) and (4) 
quantified in Table 5 provides most probably a lower bound for the statistical 
significance of the hypothesis that both C and D should be chosen non-zero, 
above the 99.9999% confidence level. Indeed, a non-Gaussian noise with a fat- 
tailed distribution would be expected to decrease the relevance of competing 
formulas, whose performance could be scrambled and be made fuzzy. The clear 
and strong result of the Wilk tests with assumed Gaussian noises thus confirm 
a very strong significance of both formulas (1) and (4). 



2.6 Determination of tc 

The critical time defines the real starting time of the anti-bubbles and is 
an important parameter for quantifying the synchronization between different 
stock markets. It is thus important to investigate how robust is its determi- 
nation by our fitting procedure ^ . In this section, we discuss two markets to 
illustrate the typical situation, the French stock index as an example of a 
bearish anti-bubble and the Australia stock index as an example of a bullish 
anti-bubble. To test for the robustness of the determination of tc, we follow the 
analysis of [43] on the USA Standard and Poor index and take seven different 
values for igtart from Jun-01-2000 to Dec-01-2000 for each index. Figure 29 



^ Wc do not discuss here other approaches for the estimation of tc-, such as using 
Shank's transformation, the generaiized q anaiysis, the parametric fitting approach, 
and so on (see Ref. [43] and reference therein). 
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shows the seven best fits, one for each tstart; for the French stock index. The 
corresponding fitting parameters are hsted in Table 6. Three fits in Fig. 29 are 
shghtly different from the rest especiaUy in the early days of the anti-bubble. 
The starting dates of these three fits are Oct-01-2000, Nov-01-2000 and Dec- 
01-2000. The slightly different nature of these three cases is also reflected in 
discernable variations in the fitting parameters listed in Table 6. In particular, 
they identify a critical tc at the end of October, 2000 rather than mid- August, 
2000. They also have slightly larger exponents a, lower log- frequencies uj and 
smaller A. Despite these differences, all the fits are quite robust indicating a 
critical time at or shghtly after August, 2000. 

Figure 30 shows the seven best fits, one for each tstart, for the Australian stock 
index. The relevant parameters are listed in Table 7. We also observe three 
fits starting on Jun-01-2000, Nov-01-2000 and Dec-01-2000 that have relatively 
later predicted tc, slightly larger exponents a, lower log-frequencies cu. But the 
predicted index value at tc (i.e., e^) are almost the same. The critical time 
tc of the Australian anti-bubble is also clustered around mid- August, 2000. 
This is consistent with the fact that both including extra data earlier than tc 
(istart ^ tc) and truncating data after tc (istart ^ tc) will reduce the precision 
of the determination of tc and deteriorate the quality of fits. 

We nevertheless have to note that not all the indices give such robust re- 
sults. The log-periodic oscillations in the initial days of some anti-bubbles are 
completely spoiled by noise, where different effects overwhelm the herding be- 
havior thought to be at the origin of the log-periodic power law patterns. For 
instance, the existence of an anti-bubble in the Peruvian stock index is quite 
questionable in view of the particularities in its fitting parameters. 



3 Correlation across different markets and synchronization of the 
anti-bubbles 

One of the most remarkable results obtained so far is that most of the anti- 
bubbles started between August and November, 2000, with very similar time 
evolutions as quantified by the formulas (1) and (4) and by the Tables 1 and 
2. This suggests that the triggering of almost simultaneously occurring anti- 
bubbles is an international event. Figure 31 summarizes our main message by 
superimposing the stock market indices of seven countries (US S&P 500, the 
Netherlands, Prance, Germany, Norway, UK, Spain). The ordinate plots the 
normalized values \p{t) — (p)]/(Jp of each index as a function of time, where 
{p) is the mean whose substraction accounts for a country-specific translation 
in price and cxp is the standard deviation for each index which accounts for 
a country-specific adjustment of scale. This remarkable collapse onto a single 
master curve does not rely on any parametric fit. It demonstrates maybe more 
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clearly than by any other means the extraordinary strong synchronization of 
the anti-bubble regime in the major western markets. Other markets exhibit 
a higher variability and have not be represented on this curve for clarity. 

It is well-known that the October 1987 crash was an international event, oc- 
curring within a few days in all major stock markets [1]. It is also often been 
noted that smaller West-European stock markets as well as other markets 
around the world are influenced by dominating trends on the USA market. 
In this spirit, in [17], a set of secondary stock markets were shown to ex- 
hibit well-correlated "anti-bubbles" triggered by a rash of crises on emerging 
markets in early 1994. In this case, the synchronization occurred between 
West-European markets which were decoupled from the USA markets. This 
suggests that smaller stock markets can weakly synchronize not only because 
of the over-arching influence of the USA market, but also independent of the 
USA market due to external factors such as the Asian crisis of 1994. 

Here, we have shown the occurrence of the synchronization of a large majority 
of markets with significant volumes into a collective anti-bubble, that includes 
the USA markets, most of the European markets as well as the developed 
Asian markets and a few other markets worldwide (see the list given in table 
1). Motivated by this result, we turn now to a series of non-parametric tests 
exploring the nature and amplitude of this worldwide synchronization, in order 
to attempt to cast addition light on this remarkable event. 

In the following, we investigate several measures of correlation, or more gen- 
erally of inter-dependence, between each index and the USA S&P500 index 
taken as a reference, in order to test whether we could have otherwise detected 
the synchronization unravelled by our log-periodic analysis. These measures of 
inter-dependence use the cross-correlation of weekly returns, linear regressions 
of indices and of their returns, a synchronization ratio of joint occurrences of 
ups and downs and an event synchronization method recently introduced [27]. 
These different measures confirm that the inter-dependence between the ma- 
jor western markets has slightly increased as a function of time in the last 
decade and especially since the Fall of 2000, confirming weakly the qualitative 
message contained in our results of the occurrence of a synchronized anti- 
bubble worldwide. However, these more standard measures of dependence do 
not come near the log-periodic analysis in the strength of the signal. 



3.1 Cross- correlation of weekly returns 

The formulas (1) and (4) have been applied to the prices and the insight into 
the existence of a synchronization comes from a comparison between these fits 
on the index prices. In order to study the cross-correlation between different 
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indices, we need to study the index returns which are approximately stationary 
thus ensuring reasonable convergence properties of the correlation estimators. 
We thus follow a procedure similar to that of Ref. [30] for the estimation of 
the cross-correlation coefficients of monthly percentage changes in major stock 
market indexes from June 1981 to September 1987. 

In order to minimize noise, we smooth the price time series with a causal 
Savitzky-Golay filter with eight points to the left of each point ("present 
time"), zero point to its right and a fourth order polynomial [26]. This provides 
a smoothed price time series p{t). We then construct the return time series 
and then obtain the cross-correlation functions of p{t). We use weekly returns, 
as a compromise between daily and monthly returns to minimize noise and 
maximize the data set size. The weekly returns are defined on the smoothed 
price time series p(t) as 

r{t) = ]n\p{t)/Pit - 7)] . (9) 

We calculate the correlation coefficients C (t) of the stock indices in a moving 
window of 65 trading days (or about a quarter in calendar days). We present 
our results obtained for the cross- correlation between the USA S&P500 index 
and the stock market indices of the Netherlands, Prance, Japan, Germany, 
UK, Hong Kong, Australia, Russia and China, which are typical. 

As illustrated in Fig. 32, the European markets have rather strong correlations 
with the American market with an average correlation coefficient of 0.44±0.05. 
The cross-correlation coefficients of the smoothed weekly returns for Hong 
Kong, Australia, Russia and China are shown in Fig. 33 as a function of 
time. Their average cross-correlation coefficients are relatively weaker than 
those for the European markets, with values respectively equal to 0.30 ±0.05, 
0.34 ± 0.05, 0.35 ± 0.05 and 0.21 ± 0.05 for Japan, Hong Kong, Austraha and 
Russia. The average cross-correlation coefficient for China is shghtly negative 
(—0.06 ± 0.05), indicating that the Chinese stock market seems practically 
uncorrclated from the western markets. The uncertainties and fluctuations of 
the variables C{t) are determined by a bootstrap simulation of 1000 series of 
reshuffled returns which gives a standard deviation ac = 0.12. 

Interestingly, Fig. 32 shows that the cross- correlation coefficients of the Eu- 
ropean markets with the American market increases slowly with time. This 
property is weaker for Japan and Hong Kong and is completely absent for 
Russia and China. While qualitatively compatible, the evidence for a slow 
increase of the cross-correlation is not sufficiently precise to relate precisely 
to our previous finding of a strong synchronization of an anti-bubble regime 
since the summer of 2000. 

To refine the evidence for an increase in correlation, we investigate the cor- 
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relation between the USA S&P500 index and nine other indices (Netherlands 
(HL), France (FR), Japan (JP), Germany (DE), United Kingdom (UK), Hong 
Kong (HK), Australia (AU), Russia (RU), and China (CN)) in two periods, 
[Jun-04-1997, Aug-09-2000] and [Aug-10-2000, Sep-04-2002]. Table 8 shows 
the P coefficients and corresponding correlation coefficients 7 of the weekly 
returns in these two periods. The coefficients /3i and /32 for the two periods of 
each index are obtained by using the well-known linear regression of the time 
series of returns of each index against the time series of returns of the S&P500 
index. Such an approach led Roll [30] to conclude on the existence of a partic- 
ularly strong synchronization during and after the crash of Oct. 1987 seen by 
the fact that the beta's of the different indices against a world market index 
were anomalously large. The two correlation coefficients 71 and 72 are directly 
evaluated. The values of the slope (3 [30] and the linear correlation coefficient 
7 are hsted in Table 8. Fig. 34 plots the returns of four European indices as 
a function of the returns of the S&P500 for each of the two periods. This 
figure and Table 8 confirm a significant increase of the correlations from the 
period [Jun-04-1997, Aug-09-2000] to the period [Aug-10-2000, Sep-04-2002]. 
Fig. 35 plots the returns of the indices of Hong Kong, Australia, Russia and 
China as a function of the returns of the S&P500 for each of the two periods. 
Table 8 and Fig. 35 show a significant increase in correlation from the period 
[Jun-04-1997, Aug-09-2000] to the period [Aug-10-2000, Sep-04-2002] only for 
Hong Kong and Australia. Russia gives a marginal signal and China none. 

Table 8 and Fig. 34 clearly confirm a strong increase in the correlation between 
the USA stock market and the European indices and some non-European in- 
dices from the period [Jun-04-1997, Aug-09-2000] to the period [Aug-10-2000, 
Sep-04-2002] , in agreement with the evidence of the log-periodic synchroniza- 
tion documented above. 



3.2 Synchronization of weekly returns 

We now discuss another intuitive measure for the characterization of the syn- 
chronization of weekly returns between different world stock markets. We use 
a moving window, whose size is 65 trading days, corresponding to 13 weeks. In 
this moving window, we define the synchronization factor R(t) as the fraction 
of weeks among the 13 weeks for which a given index return has the same 
sign as that of the S&P500 index. By definition, < R{t) < 1. R{t) = 1 
(respectively R{t) = 0) corresponds to full synchronization (respectively per- 
fect anti-synchronization). R{t) — 0.5 corresponds to independent time series 
whose weekly returns have mutually random signs. 

We calculate the synchronization factor R{t) between the USA S&P500 in- 
dex and the indices of the Netherlands, France, Japan, Germany, UK, Hong 
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Kong, Australia, Russia and China. As shown in Figs. 36 and 37, all con- 
sidered indices have R{t) significantly larger than 0.5 except for China for 
which R{t) = 0.49. The uncertainties and fluctuations of the variables R{t) 
are determined by a bootstrap simulation of 1000 series of reshuffled returns 
which give a standard deviation ctr = 0.06. Again, the European markets have 
relatively higher synchronization factors and their R{t) increase clearly with 
time. Not only R{t) is consistently at its highest long-term average level in the 
last few years for all markets, except for Russia and China, we can also note a 
very strong and significant increase of R{t) over the last year with much less 
fiuctuations. Only Russia and China among the eight indices escape from this 
world-wide synchronization. 



3.3 Time resolved event synchronization of the index time series 



In view of the importance of characterizing the dependence between different 
markets, we present yet another measure of the synchronization of weekly 
returns across different world stock markets. The two previous analyses were 
based on the time series of weekly returns. The present analysis measures the 
synchronization between different index time series by quantifying the relative 
timings of specific events in the time series, following the algorithm initially 
introduced in [27]. 

Given two index time series ^\t) and we define "events" as large 

market velocities. We define the market velocity as a coarse-grained measure 
of the slope of the price as a function of time. To obtain this coarse-grained 
measure, we apply on the prices p{t) a causal Savitzky-Golay fourth-order 
polynomial filter with eight points on the left and no point on the right. 
The velocity v{t) at time t is define as the analytical time derivative of the 
coarse-grained p{t). "Large velocities" are defined by the condition 1^(^)1 > d, 
where d is a threshold chosen here equal to d — 0.001. The times when the 
velocities v^^\t) and v^'^\t) of the two index time series p^^\t) and p^'^\t) obey 
the condition \v(t)\ > d are denoted respectively {i = 1, . . . ,mi) and t^^^ 
(j = 1, . . . , m2). The degree of synchronization is then quantified by counting 
the number of times an event (l'y(t)l > d) appears in time series p^^^t) shortly 
after it appears in time series p*^^' (t) . This number is estimated by the following 
formula 

nil ni2 

c(i|2) = EE >^^. (10) 
i=i j=i 
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with 



1 if < - tf < Tij 



1/2 if = , (11) 



otherwise 



where 

= min {e - t?, - t% - tf \ tf - t^^ . (12) 



Likewise, c(2|l) is calculated in a similar manner. The symmetrical combina- 
tion 

C = !lM±£(3i), (13) 

mim2 



called the synchronization index, measures the synchronization of the events 
and thus of the two time series. By construction, < Q < 1- The cases of 
Q = 1 and Q = correspond respectively to full synchronization and absence 
of synchronization of events of the two index time series. 

We calculate Q{t) in a moving window of 65 trading days as before between 
the USA S&P500 index on the one hand and the stock markets of the Nether- 
lands, France, Japan, Germany, UK, Hong Kong, Australia, Russia and China 
on the other hand. Fig. 38 and Fig. 39 shown that all stock markets have Q{t) 
significantly larger than 0.5 except for China. Again, an increasing trend ap- 
pears clearly for the European markets. Furthermore, the period since the 
winter of 2000 has significantly larger Q{t) compared with the earlier time for 
the European markets, Hong Kong and Australia. Q{t) is especially large and 
regular for The Netherland (HL) since mid-1999 and all other markets also 
have a very high synchronization index Q{t) since 2001. 



4 Discussion 



Following our previous investigation of the USA Standard and Poor index anti- 
bubble that started in August 2000 [43], we have analyzed the major stock 
market indices worldwide and found that a vast majority of European and 
Western countries as well as many other indices exhibit practically the same 
log-periodic power law anti-bubble structure as found for the USA S&P500 
index. In addition, these anti-bubbles are found to start approximately at the 
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same time, August 2000, in all these markets. This shows a remarkable degree 
of synchronization worldwide which, to our knowledge, has never been seen at 
any other time other than during and in the (short) aftermath of the October 
1987 crash. 

To test further this synchronization, we have also used several standard and 
less standard measures of correlation, dependence and synchronization be- 
tween the USA S&P500 index and other world markets. These measures con- 
firm the existence of significant increase of dependence in the last decade and 
still a larger increase in the last one-two years. However, these measures come 
nowhere close to the clarity of the signal of the extraordinary strong synchro- 
nization found using the log-periodic power law analysis. This is due to the 
fact that the log-periodic power law analysis is not sensitive to detailed phases 
in the oscillations (translated in shghtly shifted effective time units in different 
markets) and detects only the robust universal unit-independent discrete scale 
invariant features of the price trajectories. 

What triggered the worldwide anti-bubble in August 2000? The international 
descent of many of the worldwide stock markets since 2000 suggests the 
strengthening of globahzation and the leading impact of the USA. In this 
respect, if history is any guide, the historical record on financial crises shows 
that they are often accompanying surges of globalization in the past, includ- 
ing events as far back as in the 19th century such as during the gold standard 
period of 1880-1913 [24]. Bordo and Murshid [4] compared various charac- 
teristics of the cross-country transmission of shocks in the financial markets 
of both advanced and emerging countries during two periods of globalization 
- the pre- World War I classical gold standard era, 1880-1914, and the post- 
Bretton Woods era, 1975-2000. They found that financial market shocks were 
more globalized before 1914 compared to the present and interpret this re- 
sult by the growing financial maturity of advanced countries and the widening 
of the center to include a more diverse group of countries spanning several 
regions. Our findings temper Bordo and Murshid's results and suggest a pos- 
sible transition to a stronger integration and globalization fostered by several 
factors, including corporate and financial globalization, and the rapid develop- 
ment, adoption and use of information and communications technology. Our 
results also confirm those of Goetzmann et al. [9] who find that the correlation 
structure of the major world equity markets over 150 years vary considerably 
through time and are highest during periods of economic and financial integra- 
tion such as the late 19th and 20th centuries. Goetzmann et al. [9] also stress 
that such increase of correlation implies that diversification benefits to global 
investing relies increasingly on investment in emerging markets, in agreement 
with our results on the weaker synchronization of emerging markets. Our re- 
sults can also be seen to add to the literature on contagion, usually defined as 
correlation between markets in excess of what would be implied by economic 
fundamentals, by providing a new technical tool. 
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Table 1 

Parameters of the fits of the indices indicated in the first column using the first-order 
formula (1) from igtart to September, 30, 2002. All these indices are in the so-called 
bearish anti-bubble regime, qualified by the fact that the coefficient B is negative. 
The exponents a of the leading power law in formula (1) are found either larger or 
smaller than 1, corresponding to an accelerating (respectively decelerating) decrease 
of the prices as a function of time (see text), x denotes the root-mean-square (r.m.s.). 
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3.96 


Spain 


00/09/14 


00/10/04 


0.93 


7.52 


3. 


,15 


6.90 


-0.82 


0.31 


3.98 


Switzerland 


00/08/23 


00/11/18 


1.00 


6.76 


5. 


,16 


9.01 


-0.70 


-0.26 


3.67 


UK 


00/09/04 


00/10/23 


1.00 


7.58 


0. 


,00 


8.77 


-0.55 


-0.17 


3.17 


Israel 


00/08/28 


00/09/09 


0.18 


11.45 


4. 


,02 


6.61 


-205 


17.7 


4.51 


Brazil 


00/08/14 


00/08/12 


0.87 


10.15 


4. 


,74 


9.74 


-1.36 


0.56 


6.20 


Hong Kong 


00/07/21 


00/02/26 


0.90 


7.65 


1. 


,95 


9.90 


-1.69 


-0.29 


5.27 


India 


00/07/12 


00/05/26 


0.78 


6.34 


5. 


,61 


8.47 


-2.92 


-0.85 


4.89 


Peru 


00/09/11 


99/10/28 


0.15 


19.80 


1. 


,16 


7.45 


-118 


-20.9 


3.22 


Taiwan 


00/07/17 


00/01/24 


0.38 


9.22 


1. 


,01 


9.26 


-63.0 


17.7 


8.50 


Czech 


00/07/28 


00/08/13 


0.53 


4.61 


1. 


,19 


6.32 


-12.5 


5.20 


4.80 


Argentina 


00/07/14 


00/05/02 


1.54 


7.25 


5. 


,28 


6.23 


-0.03 


-0.02 


12.2 


Turkey 


00/07/10 


00/08/28 


0.13 


9.56 


2. 


,88 


9.62 


-160 


63.7 


11.8 
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Table 2 

Values of the parameters of the fits to the 21 indices of table 1 given in the first 
column using formula (4) from fgtart to September, 30, 2002. As in table 1, the 
bearish nature of the markets is identified by the fact that B < 0. Note that the 
amplitude D of the harmonics is often very significant compared with the amplitude 
C of the mean angular log-frequency. 



Stock 


/ sla.rl 


/. 


c\ 


jj 


Ol 


02 


A 


lif'D 


10'^ c 


UfD 


lO^x 


Netherlands 


00/09/04 


00/07/20 


0.80 


11.07 


4.75 


0.91 


6.60 


-2.82 


0.64 


0.32 


3.72 


France 


00/09/04 


00/07/04 


0.71 


11.24 


0.28 


4.25 


8.90 


-5.76 


-0.91 


0.48 


3.47 


USA DowJoncs 


00/09/06 


00/06/18 


0.72 


12.09 


0.83 


2.16 


9.33 


-1.45 


-0.65 


-0.34 


3.00 


USA NASDAQ 


00/08/20 


00/09/02 


0.26 


10.00 


3.37 


0.00 


8.77 


-269.2 


-24.73 


-9.18 


5.70 


Japan 


00/08/28 


00/08/06 


0.77 


7.74 


3.35 


0.00 


9.74 


-3.72 


-0.89 


0.21 


3.69 


Belgium 


00/11/06 


00/06/21 


1.14 


12.24 


2.67 


0.49 


8.05 


-0.12 


0.05 


-0.02 


2.83 


Denmark 


00/10/24 


00/05/04 


0.64 


13.35 


0.87 


0.00 


6.04 


-7.51 


1.05 


0.36 


2.86 


Germany 


00/09/04 


00/10/06 


0.94 


8.47 


3.61 


4.58 


8.86 


-1.20 


0.41 


0.12 


3.96 


Norway 


00/09/05 


00/07/13 


0.87 


10.92 


2.03 


1.78 


6.82 


-1.61 


-0.45 


0.13 


3.58 


Spain 


00/09/14 


99/07/03 


0.88 


15.70 


1.75 


1.28 


7.14 


-1.14 


-0.18 


0.08 


3.46 


Switzerland 


00/08/23 


00/11/23 


0.95 


6.72 


5.54 


0.00 


9.01 


-0.98 


-0.35 


-0. 07 


3.52 


UK 


00/09/04 


00/07/21 


0.84 


10.78 


3.56 


1.52 


8.83 


-1.53 


-0.30 


-0.16 


2.82 


Israel 


00/08/28 


00/09/09 


0.19 


11.41 


4.26 


0.00 


6.60 


-194. 


17.28 


-1.73 


4.49 


Brazil 


00/08/14 


00/02/23 


1.08 


7.99 


6.24 


4.69 


9.83 


-0.34 


-0.07 


-0.08 


5.41 


HongKong 


00/07/21 


00/01/30 


0.41 


7.52 


2.21 


0.41 


10.39 


-71.68 


-4.88 


3.03 


4.63 


India 


00/07/12 


99/12/04 


0.03 


8.84 


5.96 


0.73 


17.34 


-7544 


-63 


39.5 


4.12 


Peru 


00/09/11 


00/02/26 


0.54 


7.69 


2.98 


4.99 


7.19 


-1.67 


1.20 


1.92 


2.76 


Taiwan 


00/07/17 


99/05/07 


0.31 


7.34 


1.93 


4.44 


9.95 


-181.0 


-14.51 


21.50 


5.41 


Czech 


00/07/28 


00/08/11 


0.40 


4.59 


1.27 


0.77 


6.39 


-33.90 


9.69 


-4.46 


3.94 


Argentina 


00/07/14 


99/08/23 


2.00 


11.23 


0.67 


4.40 


6.25 


-0.00 


-0.00 


0.00 


8.20 


Turkey 


00/07/10 


00/08/28 


0.20 


9.55 


2.94 


1.32 


9.52 


-78.16 


47.12 


18.47 


11.0 
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Table 3 

Values of the parameters of the fits of the indices indicated in the first column using 
the first-order formula (1) from fgtart to September, 30, 2002. The bullish nature 
of the markets refers to the the fact that B > 0. The exponents a of the leading 
power law take values cither greater or less than 1 indicating either an upward or a 
downward overall large scale curvature of the price trajectory. 



Stock 


^stcirt 


tc 


a 


u 




A 


10^5 


lO^C 




Australia 


00/09/04 


00/08/08 


0.77 


10.85 


3.55 


8.08 


0.07 


-0.37 


1.76 


Mexico 


00/09/04 


00/08/08 


0.80 


10.11 


4.42 


8.70 


0.32 


0.86 


4.02 


Indonesia 


00/09/16 


00/08/09 


1.06 


10.30 


5.83 


5.99 


0.09 


-0.21 


3.90 


Korea 


00/10/12 


00/11/15 


1.37 


6.69 


5.67 


6.29 


0.05 


-0.05 


5.00 


Thailand 


00/09/20 


00/08/16 


0.93 


9.45 


5.16 


5.61 


0.58 


-0.39 


4.37 


Russia 


00/12/01 


00/10/08 


0.92 


7.93 


3.33 


7.24 


2.76 


-0.68 


4.85 



Table 4 

Values of the parameters of the fits of the indices indicated in the first column 
using formula (4) from tgtart to September, 30, 2002. 



Stock 


^stcirt 


tc 


a 


UJ 


4>i 


<l>2 


A 


lO^S 


lO^C 






Australia 


00/09/04 


00/08/10 


0.64 


10.89 


3.42 


0.51 


8.08 


0.26 


-0.80 


-0.26 


1.53 


Mexico 


00/09/04 


00/06/27 


0.73 


5.97 


4.23 


1.65 


8.69 


0.72 


0.39 


-1.15 


3.91 


Indonesia 


00/09/16 


00/11/16 


0.96 


3.47 


1.33 


3.10 


6.04 


-0.25 


-0.50 


-0.62 


2.71 


Korea 


00/10/12 


00/08/10 


1.00 


5.04 


1.90 


1.90 


6.24 


0.46 


-0.23 


-0.25 


4.28 


Thailand 


00/09/20 


00/07/27 


0.43 


5.18 


3.62 


2.22 


5.48 


20.83 


3.90 


6.62 


3.54 


Eussia 


00/12/01 


00/10/12 


0.91 


7.85 


3.90 


0.17 


7.23 


3.07 


-0.75 


0.10 


4.57 
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Table 5 

Likelihood-ratio (Wilk) test of hypothesis Hi against Hq and of hypothesis H2 
against Hi, where Hq is the hypothesis that C = in (1) (pure power law fit), Hi 
is the hypothesis that D = in (4) (corresponding to the log-periodic function (1) 
without any harmonics) and H2 is the hypothesis that D ^ 0, corresponding to the 
log-periodic function (4) which includes an harmonic at 2io. Each stock index time 
series whose country/area is given in the first column starts at igtart and ends at 
September 30, 2002. The column n gives the number of the data points for each fit. 
The cjj's are the standard deviations for hypothesis Hj with j = 0, 1, 2 of the fits to 
the data associated with formula (1) and (4). The Tjj^i are the log-likelihood-ratios 
defined in expression (8). The Pjj+i are the probability for exceeding Tjj^i under 
the assumption of that hypothesis Hj holds. 



Stock 


^staxt 


n 




0-1 


0"2 


To,i 


Tl,2 


Po,i(%) 


Pl,2(%) 


Netherlands 


00/09/04 


508 


0.078 


0.043 


0.037 


609.2 


148.5 


< 


10- 


-4 


< 10- 


-4 


Prance 


00/09/04 


520 


0.070 


0.040 


0.035 


589.0 


144.2 


< 


10- 


-4 


< 10- 


-4 


USA DowJones 


00/09/06 


500 


0.053 


0.034 


0.030 


459.8 


118.1 


< 


10- 


-4 


< 10- 


-4 


USA NASDAQ 


00/08/20 


511 


0.099 


0.064 


0.057 


456.2 


111.5 


< 


10- 


-4 


< 10- 


-4 


Japan 


00/08/28 


498 


0.069 


0.040 


0.037 


551.7 


70.7 


< 


10- 


-4 


< 10- 


-4 


Belgium 


00/11/06 


459 


0.056 


0.034 


0.028 


466.7 


158.6 


< 


10^ 


-4 


< 10- 


-4 


Denmark 


00/10/24 


456 


0.052 


0.032 


0.029 


457.1 


92.2 


< 


10- 


~4 


< 10- 


-4 


Germany 


00/09/04 


504 


0.083 


0.044 


0.040 


637.9 


105.7 


< 


10- 


-4 


< 10- 


-4 


Norway 


00/09/05 


543 


0.081 


0.040 


0.036 


772.0 


110.8 


< 


10- 


-4 


< 10- 


-4 


Spain 


00/09/14 


483 


0.062 


0.040 


0.035 


421.9 


133.0 


< 


10- 


-4 


< 10- 


-4 


Switzerland 


00/08/23 


505 


0.062 


0.037 


0.035 


526.9 


40.9 


< 


10- 


-4 


< 10- 


-4 


UK 


00/09/04 


506 


0.050 


0.032 


0.028 


464.5 


120.4 


< 


10- 


-4 


< 10- 


-4 


Israel 


00/08/28 


390 


0.058 


0.045 


0.045 


194.0 


2.7 


< 


10- 


-4 


10.03 


Brazil 


00/08/14 


505 


0.094 


0.062 


0.054 


418.9 


138.0 


< 


10" 


-4 


< 10- 


-4 


HongKong 


00/07/21 


523 


0.070 


0.053 


0.046 


298.3 


135.6 


< 


10- 


-4 


< 10- 


-4 


India 


00/07/12 


528 


0.073 


0.049 


0.041 


428.7 


180.2 


< 


10- 


-4 


< 10- 


-4 


Peru 


00/09/11 


488 


0.051 


0.032 


0.028 


438.3 


151.9 


< 


10- 


-4 


< 10- 


-4 


Taiwan 


00/07/17 


518 


0.114 


0.085 


0.054 


307.3 


466.7 


< 


10- 


-4 


< 10- 


-4 


Czech 


00/07/28 


507 


0.052 


0.048 


0.039 


72.4 


199.0 


< 


10- 


-4 


< 10- 


-4 


Argentina 


00/07/14 


502 


0.148 


0.122 


0.082 


194.0 


401.1 


< 


10- 


-4 


< 10- 


-4 


Turkey 


00/07/10 


534 


0.146 


0.118 


0.110 


233.1 


67.8 


< 


10- 


-4 


< 10- 


-4 


Australia 


00/09/04 


509 


0.031 


0.018 


0.015 


588.4 


144.1 


< 


10- 


-4 


< 10- 


-4 


Mexico 


00/09/04 


493 


0.078 


0.040 


0.039 


655.4 


25.9 


< 


10- 


-4 


« 10- 


-4 


Indonesia 


00/09/16 


454 


0.080 

26 
0.108 


0.039 


0.027 


658.6 


329.1 


< 


10- 


-4 


< 10- 


-4 


Korea 


00/10/12 


456 


0.050 


0.043 


705.2 


142.3 


< 


10- 


-4 


< 10- 


-4 


Thailand 


00/09/20 


481 


0.076 


0.044 


0.035 


535.1 


203.0 


< 


10- 


-4 


< 10- 


-4 



Table 6 

Fitting parameters with equation (1) for the French market stock index with dif- 
ferent tstait indicated in the first column. 



^stcirt 


tc 


a 


UJ 




A 


B 


c 


X 


00/06/01 


00/08/20 


0.91 


9.23 


4.26 


8.81 


-0.00154 


0.00035 


0.0387 


00/07/01 


00/08/15 


0.89 


9.53 


5.48 


8.81 


-0.00168 


-0.00038 


0.0389 


00/08/01 


00/08/19 


0.89 


9.36 


3.45 


8.81 


-0.00165 


0.00038 


0.0393 


00/09/01 


00/08/30 


0.91 


8.90 


3.37 


8.79 


-0.00145 


-0.00035 


0.0399 


00/10/01 


00/10/23 


0.97 


7.56 


3.19 


8.74 


-0.00101 


0.00029 


0.0400 


00/11/01 


00/10/25 


0.98 


7.46 


0.70 


8.74 


-0.00096 


-0.00027 


0.0407 


00/12/01 


00/10/23 


1.05 


7.48 


0.50 


8.72 


-0.00061 


-0.00018 


0.0413 



Table 7 

Fitting parameters with equation (1) for the Australian market stock index with 
different tstart indicated in the first column. 



^stcirt 


tc 


a 


CO 


4> 


A 


B 


C 


X 


00/06/01 


00/09/23 


0.76 


9.41 


0.88 


8.08 


0.00002 


-0.00044 


0.0173 


00/07/01 


00/08/21 


0.80 


10.25 


4.42 


8.09 


0.00001 


0.00032 


0.0173 


00/08/01 


00/08/20 


0.78 


10.29 


4.17 


8.09 


0.00002 


0.00036 


0.0176 


00/09/01 


00/08/14 


0.77 


10.57 


2.29 


8.08 


0.00006 


0.00037 


0.0176 


00/10/01 


00/08/11 


0.78 


10.68 


4.67 


8.08 


0.00007 


-0.00035 


0.0177 


00/11/01 


00/09/26 


0.72 


9.35 


4.47 


8.08 


0.00003 


0.00055 


0.0180 


00/12/01 


00/09/11 


0.67 


9.78 


4.66 


8.08 


0.00007 


-0.00073 


0.0181 
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Table 8 

/3's and correlation coefficients 7 between the USA S&P500 index and nine indices 
(Netherlands (HL), France (FR), Japan (JP), Germany (DE), United Kingdom 
(UK), Hong Kong (HK), Australia (AU), Russia (RU), and China (CN)) in two 
periods. Period 1 is [Jun-04-1997, Aug-09-2000] and period 2 is [Aug-10-2000, Sep- 
04-2002] . The return is the logarithm of the ratio between two successive prices with 
time lag of 30 trading days. 



INDEX 


HL 


FR 


JP 


DE 


UK 


HK 


AU 


RU 


CN 




1.07 


0.93 


0.98 


0.92 


0.68 


1.34 


0.60 


2.21 


-0.05 


71 


0.73 


0.67 


0.62 


0.83 


0.49 


0.51 


0.68 


0.42 


0.03 




1.06 


1.03 


1.17 


0.76 


0.66 


0.92 


0.45 


0.99 


0.02 


72 


0.87 


0.90 


0.86 


0.87 


0.59 


0.72 


0.80 


0.46 


0.02 
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Fig. 1. The stock market index of Netherlands (fine noisy line) and its fits with the 
simple log-periodic formula (1) (thick line) and with formula (4) incorporating an 
harmonics at 2u! (dashed line). The starting date for the fits is tstart = 04-Sep-2000 
and the ending date of the fit is September, 30, 2002. The parameter values of 
the fit with (1) are tc = 28-Aug-2000, a = 1.05, oj = 9.16, (j) = 1.63, A = 6.53, 
B = -0.00055, and C = -0.00018. The r.m.s. of the fit errors is 0.043. The param- 
eter values of the fit with (4) are tc = 20-Jul-2000, a = 0.80, oj = 11.07, = 4.75, 
(/>2 = 0.91, A = 6.60, B = -0.00282, C = 0.00064, and D = 0.00032. The r.m.s. of 
the fit errors is 0.0372. The formula including the harmonics reduces the r.m.s. of 
fit errors by 13.6%. 
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Fig. 2. The stock market index of France (fine noisy line) and its fits with the 
simple log-periodic formula (1) (thick line) and with formula (4) incorporating an 
harmonics (dashed line). The starting date for the fits is igtart = 04-Sep-2000. The 
parameter values of the fit with (1) are tc = 30-Aug-2000, a = 0.92, uj = 8.88, 
(/> = 3.54, A = 8.79, B = -0.00139, and C = -0.00033. The r.m.s. of the fit errors 
is 0.0399. The parameter values of the fit with (4) are tc = 04-Jul-2000, a = 0.71, 
uj = 11.24, 01 = 0.28, (1)2 = 4.25, A = 8.90, B = -0.00576, C = -0.00091, and 
D = 0.00048. The r.m.s. of the fit errors is 0.0347. The formula including the 
harmonics reduces the r.m.s. of fit errors by 12.9%. 
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Fig. 3. The stock market DowJones index of USA (fine noisy line) and its fits with 
the simple log-periodic formula (1) (thick line) and with formula (4) incorporating 
an harmonics (dashed line). The starting date for the fits is tstart = 06-Sep-2000. 
The parameter values of the fit with (1) are tc = 15-Aug-2000, a = 1.05, lo = 9.76, 
(p = 4.01, A = 9.30, B = -0.00017, and C = -0.00011. The r.m.s. of the fit 
errors is 0.0337. The parameter values of the fit with (4) are tc = 18-Jun-2000, 
a = 0.72, uj = 12.09, = 0.83, 02 = 2.16, A = 9.33, B = -0.00145, C = -0.00065, 
and D = —0.00034. The r.m.s. of the fit errors is 0.03. The formula including the 
harmonics reduces the r.m.s. of fit errors by 11.1%. 
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Fig. 4. The stock market NASDAQ index of USA (fine noisy line) and its fits with 
the simple log-periodic formula (f) (thick line) and with formula (4) incorporating 
an harmonics (dashed line). The starting date for the fits is tstart = 20-Aug-2000. 
The parameter values of the fit with (1) are = 02-Sep-2000, a = 0.26, uj = 10.00, 
(j) = 3.37, A = 8.74, B = -0.25128, and C = -0.02360. The r.m.s. of the fit errors 
is 0.0636. The parameter values of the fit with (4) are tc = 02-Sep-2000, a = 0.26, 
UJ = 10.00, 01 = 3.37, </>2 = 0.00, A = 8.77, B = -0.26922, C = -0.02473, and 
D = —0.00918. The r.m.s. of the fit errors is 0.057. The formula including the 
harmonics reduces the r.m.s. of fit errors by 10.3%. 
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Fig. 5. The stock market index of Japan (fine noisy line) and its fits witli the 
simple log-periodic formula (1) (thick line) and with formula (4) incorporating an 
harmonics (dashed line). The starting date for the fits is tstart = 28-Aug-2000. The 
parameter values of the fit with (1) are tc = 06-Aug-2000, a = 0.79, u = 7.74, 
4> = 3.34, A = 9.74, B = -0.00340, and C = -0.00086. The r.m.s. of the fit 
errors is 0.0396. The parameter values of the fit with (4) are tc = 06-Aug-2000, 
a = 0.77, uj = 7.74, (pi = 3.35, (j)2 = 0.00, A = 9.74, B = -0.00372, C = -0.00089, 
and D = 0.00021. The r.m.s. of the fit errors is 0.0369. The formula including the 
harmonics reduces the r.m.s. of fit errors by 6.9%. 
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Fig. 6. The stock market index of Belgium (fine noisy line) and its fits with the 
simple log-periodic formula (1) (thick line) and with formula (4) incorporating an 
harmonics (dashed line). The starting date for the fits is tstart = 06-Nov-2000. The 
parameter values of the fit with (1) are tc = 25-Jun-2000, a = 1.52, uj = 12.20, 
(j) = 2.85, A = 8.02, B = -0.00001, and C = 0.00000. The r.m.s. of the fit errors 
is 0.0336. The parameter values of the fit with (4) are tc = 21-Jun-2000, a = 1.14, 
UJ = 12.24, (t>i = 2.67, (1)2 = 0.49, A = 8.05, B = -0.00012, C = 0.00005, and 
D = —0.00002. The r.m.s. of the fit errors is 0.0283. The formula including the 
harmonics reduces the r.m.s. of fit errors by 15.9%. 
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Fig. 7. The stock market index of Denmark (fine noisy line) and its fits with the 
simple log-periodic formula (1) (thick line) and with formula (4) incorporating an 
harmonics (dashed line). The starting date for the fits is istart = 24-Oct-2000. The 
parameter values of the fit with (1) are tc = 03-May-2000, a = 0.78, uj = 13.37, 
4> = 0.62, A = 5.98, B = -0.00273, and C = 0.00047. The r.m.s. of the fit errors 
is 0.0316. The parameter values of the fit with (4) are tc = 04-May-2000, a = 0.64, 
uj = 13.35, = 0.87, (t>2 = 0.00, A = 6.04, B = -0.00751, C = 0.00105, and 
D = 0.00036. The r.m.s. of the fit errors is 0.0286. The formula including the 
harmonics reduces the r.m.s. of fit errors by 9.6%. 
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Fig. 8. The stock market index of Germany (fine line) and its fits witli the simple 
log-periodic formula (f ) (thick line) and with formula (4) incorporating an harmon- 
ics (dashed line). The starting date for the fits is tstart = 04-Sep-2000. The param- 
eter values of the fit with (1) are tc = 31-Aug-2000, a = 1.05, uj = 9.02, 4> = 5.66, 
A = 8.87, B = -0.00057, and C = 0.00019. The r.m.s. of the fit errors is 0.0439. 
The parameter values of the fit with (4) are tc = 06-Oct-2000, a = 0.94, uj = 8.47, 
01 = 3.61, 02 = 4.58, A = 8.86, B = -0.00120, C = 0.00041, and D = 0.00012. 
The r.m.s. of the fit errors is 0.0396. The formula including the harmonics reduces 
the r.m.s. of fit errors by 10.0%. 
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Fig. 9. The stock market index of Norway (fine line) and its fits with the simple 
log-periodic formula (1) (thick line) and with formula (4) incorporating an harmon- 
ics (dashed line). The starting date for the fits is tstart = 05-Sep-2000. The param- 
eter values of the fit with (1) are tc = 02-Oct-2000, a = 1.02, to = 8.21, (p = 4.77, 
A = 6.75, B = -0.00060, and C = 0.00022. The r.m.s. of the fit errors is 0.0396. 
The parameter values of the fit with (4) are tc = 13-Jul-2000, a = 0.87, iv = 10.92, 
01 = 2.03, 02 = 1.78, A = 6.82, B = -0.00161, C = -0.00045, and D = 0.00013. 
The r.m.s. of the fit errors is 0.0358. The formula including the harmonics reduces 
the r.m.s. of fit errors by 9.7%. 
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Fig. 10. The stock market index of Spain (fine noisy line) and its fits with the 
simple log-periodic formula (1) (thick line) and with formula (4) incorporating an 
harmonics (dashed line). The starting date for the fits is istart = 14-Sep-2000. The 
parameter values of the fit with (1) are tc = 04-Oct-2000, a = 0.93, uj = 7.52, 
(j) = 3.15, A = 6.90, B = -0.00082, and C = 0.00031. The r.m.s. of the fit errors 
is 0.0398. The parameter values of the fit with (4) are tc = 03-Jul-1999, a = 0.88, 
uj = 15.70, 01 = 1.75, (t>2 = 1.28, A = 7.14, B = -0.00114, C = -0.00018, and 
D = 0.00008. The r.m.s. of the fit errors is 0.0346. The harmonic formula reduces 
the r.m.s. of fit errors by 12.9%. 
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Fig. 11. The stock market index of Switzerland (fine noisy line) and its fits with 
the simple log-periodic formula (1) (thick line) and with formula (4) incorporating 
an harmonics (dashed line). The starting date for the fits is tstart = 23-Aug-2000. 
The parameter values of the fit with (1) are tc = 18-Nov-2000, a = 1.00, lo = 6.76, 
(j) = 5.16, A = 9.01, B = -0.00070, and C = -0.00026. The r.m.s. of the fit errors 
is 0.0367. The parameter values of the fit with (4) are tc = 23-Nov-2000, a = 0.95, 
uj = 6.72, 01 = 5.54, 4>2 = 0.00, A = 9.01, B = -0.00098, C = -0.00035, and 
D = —0.00007. The r.m.s. of the fit errors is 0.0352. The formula including the 
harmonics reduces the r.m.s. of fit errors by 4.0%. 
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Fig. 12. The stock market index of UK (fine noisy line) and its fits with the simple 
log-periodic formula (1) (thick line) and with formula (4) incorporating an har- 
monics (dashed line). The starting date for the fits is tgtart = 04-Sep-2000. The 
parameter values of the fit with (1) are tc = 23-Oct-2000, a = 1.00, u = 7.58, 
(j) = 0.00, A = 8.77, B = -0.00055, and C = -0.00017. The r.m.s. of the fit errors 
is 0.0317. The parameter values of the fit with (4) are tc = 21-Jul-2000, a = 0.84, 
uj = 10.78, (^1 = 3.56, (1)2 = 1.52, A = 8.83, B = -0.00153, C = -0.00030, and 
D = —0.00016. The r.m.s. of the fit errors is 0.0282. The formula including the 
harmonics reduces the r.m.s. of fit errors by 11.2%. 
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Fig. 13. The stock market index of Israel (fine noisy line) and its fits with the 
simple log-periodic formula (1) (thick line) and with formula (4) incorporating an 
harmonics (dashed line). The starting date for the fits is tstart = 28-Aug-2000. The 
parameter values of the fit with (1) are tc = 09-Sep-2000, a = 0.18, to = 11.45, 
(j) = 4.02, A = 6.61, B = -0.20470, and C = 0.01772. The r.m.s. of the fit errors 
is 0.0451. The parameter values of the fit with (4) are tc = 09-Sep-2000, a = 0.19, 
uj = 11.41, = 4.26, (t>2 = 0.00, A = 6.60, B = -0.19459, C = 0.01728, and 
D = —0.00173. The r.m.s. of the fit errors is 0.0449. The formula including the 
harmonics reduces the r.m.s. of fit errors by 0.4%. 
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Fig. 14. The stock market index of Brazil (fine noisy line) and its fits with the 
simple log-periodic formula (1) (thick line) and with formula (4) incorporating an 
harmonics (dashed line). The starting date for the fits is ^start = 14-Aug-2000. The 
parameter values of the fit with (1) are tc = 12-Aug-2000, a = 0.87, to = 10.15, 
(j) = 4.74, A = 9.74, B = -0.00136, and C = 0.00056. The r.m.s. of the fit errors 
is 0.062. The parameter values of the fit with (4) are tc = 23-Feb-2000, a = 1.08, 
uj = 7.99, 01 = 6.24, (t>2 = 4.69, A = 9.83, B = -0.00034, C = -0.00007, and 
D = —0.00008. The r.m.s. of the fit errors is 0.0541. The formula including the 
harmonics reduces the r.m.s. of fit errors by 12.8%. 
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Fig. 15. The stock market index of Hong-Kong (fine noisy line) and its fits with 
the simple log-periodic formula (1) (thick line) and with formula (4) incorporating 
an harmonics (dashed line). The starting date for the fits is tstart ~ 

21-Jul-2000. 

The parameter values of the fit with (1) are tc = 26-Feb-2000, a = 0.90, u = 7.65, 
(j) = 1.95, A = 9.90, B = -0.00169, and C = -0.00029. The r.m.s. of the fit 
errors is 0.0527. The parameter values of the fit with (4) are tc = 30-Jan-2000, 
a = 0.41, uj = 7.52, 4>i = 2.21, 4>2 = 0.41, A = 10.39, B = -0.07168, C = -0.00488, 
and D = 0.00303. The r.m.s. of the fit errors is 0.0463. The formula including the 
harmonics reduces the r.m.s. of fit errors by 12.2%. 
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Fig. 16. The stock market index of India (fine noisy line) and its fits with the 
simple log-periodic formula (1) (thick line) and with formula (4) incorporating an 
harmonics (dashed line). The starting date for the fits is tstart = 12-Jul-2000. The 
parameter values of the fit with (1) are tc = 26-May-2000, a = 0.78, lo = 6.34, 
(p = 5.61, A = 8.47, B = -0.00292, and C = -0.00085. The r.m.s. of the fit 
errors is 0.0489. The parameter values of the fit with (4) are tc = 04-Dec-1999, 
a = 0.03, uj = 8.84, 4>i = 5.96, 4>2 = 0.73, A = 17.34, B = -7.54442, C = -0.06297, 
and D = 0.03948. The r.m.s. of the fit errors is 0.0412. The formula including the 
harmonics reduces the r.m.s. of fit errors by 15.7%. 
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Fig. 17. The stock market index of Peru (fine noisy line) and its fits with the 
simple log-periodic formula (1) (thick line) and with formula (4) incorporating an 
harmonics (dashed line). The starting date for the fits is istart = ll-Sep-2000. The 
parameter values of the fit with (1) are tc = 28-Oct-1999, a = 0.15, u = 19.80, 
4> = 1.16, A = 7.45, B = -0.11761, and C = -0.02090. The r.m.s. of the fit 
errors is 0.0322. The parameter values of the fit with (4) are tc = 26-Feb-2000, 
a = 0.54, uj = 7.69, 4>i = 2.98, 02 = 4.99, A = 7.19, B = -0.00167, C = 0.00120, 
and D = 0.00192. The r.m.s. of the fit errors is 0.0276. The formula including the 
harmonics reduces the r.m.s. of fit errors by 14.4%. 
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Fig. 18. The stock market index of Taiwan (fine noisy line) and its fits with the 
simple log-periodic formula (1) (thick line) and with formula (4) incorporating an 
harmonics (dashed line). The starting date for the fits is tstart = 17-Jul-2000. The 
parameter values of the fit with (1) are tc = 24-Jan-2000, a = 0.38, uj = 9.22, 
(/> = 1.01, A = 9.26, B = -0.06302, and C = 0.01769. The r.m.s. of the fit errors 
is 0.085. The parameter values of the fit with (4) are tc = 07-May-1999, a = 0.31, 
UJ = 7.34, 01 = 1.93, (/)2 = 4.44, A = 9.95, B = -0.18106, C = -0.01451, and 
D = 0.02150. The r.m.s. of the fit errors is 0.0541. The formula including the 
harmonics reduces the r.m.s. of fit errors by 36.3%. 
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Fig. 19. The stock market index of Czech repubhc (fine noisy Hne) and its fits with 
the simple log-periodic formula (1) (thick line) and with formula (4) incorporating 
an harmonics (dashed line). The starting date for the fits is tstart ~ 

28-Jul-2000. 

The parameter values of the fit with (1) are tc = 13-Aug-2000, a = 0.53, lo = 4.61, 
(j) = 1.19, A = 6.32, B = -0.01246, and C = 0.00520. The r.m.s. of the fit errors 
is 0.048. The parameter values of the fit with (4) are tc = ll-Aug-2000, a = 0.40, 
uj = 4.59, = 1.27, 4>2 = 0.77, A = 6.39, B = -0.03390, C = 0.00969, and 
D = —0.00446. The r.m.s. of the fit errors is 0.0394. The formula including the 
harmonics reduces the r.m.s. of fit errors by 17.8%. 
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Fig. 20. The stock market index of Argentina (fine noisy line) and its fits with the 
simple log-periodic formula (1) (thick line) and with formula (4) incorporating an 
harmonics (dashed line). The starting date for the fits is tstart = 14-Jul-2000. The 
parameter values of the fit with (1) are tc = 02-May-2000, a = 1.54, uj = 7.25, 
(p = 5.28, A = 6.23, B = -0.00003, and C = -0.00002. The r.m.s. of the fit 
errors is 0.1223. The parameter values of the fit with (4) are tc = 23-Aug-1999, 
a = 2.00, UJ = 11.23, (/>i = 0.67, 02 = 4.40, A = 6.25, B = -0.00000, C = -0.00000, 
and D = 0.00000. The r.m.s. of the fit errors is 0.082. The formula including the 
harmonics reduces the r.m.s. of fit errors by 32.9%. 



48 



15000 




7000 1 > > > > 1 

2000 2001 2002 2003 2004 2005 

Date 

Fig. 21. The stock market index of Turkey (fine noisy line) and its fits with the 
simple log-periodic formula (1) (thick line) and with formula (4) incorporating an 
harmonics (dashed line). The starting date for the fits is tstart = lO-Jul-2000. The 
parameter values of the fit with (1) are tc = 28-Aug-2000, a = 0.13, u; = 9.56, 
(j) = 2.88, A = 9.62, B = -0.16041, and C = 0.06365. The r.m.s. of the fit errors 
is 0.1176. The parameter values of the fit with (4) are = 28-Aug-2000, a = 0.20, 
uj = 9.55, 01 = 2.94, (/)2 = 1.32, A = 9.52, B = -0.07816, C = 0.04712, and 
D = 0.01847. The r.m.s. of the fit errors is 0.1104. The formula including the 
harmonics reduces the r.m.s. of fit errors by 6.1%. 
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Fig. 22. Lomb periodograms for the 21 stock market indices shown in previous 
figures as a function of the angular log-periodic frequency u. These 21 stock market 

indices exhibit a so-called bearish anti-bubble regime characterized by i? < 0. The 
number of points in the Lomb analysis is in the range 350-530, which implies that 
the log-periodic signals are very significant for most markets with high Lomb peaks. 
The relevance of the log-periodicity is reflected in the quasi-universal value of the 
angular log-periodic frequency found in the range 7 — 10. Notice also the presence 
of an harmonic at 2lo in many of the markets outlined by the cluster of secondary 
peaks in the vicinity oi u> = 20. 
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Fig. 23. The stock market index of Australia (fine line) and its fits with the simple 
log-periodic formula (1) (thick line) and with formula (4) incorporating an har- 
monics (dashed line). The starting date for the fits is tstart = 04-Sep-2000 and the 
ending date of the fit is September, 30, 2002. The parameter values of the fit with 
(1) are tc = 08-Aug-2000, a = 0.77, uj = 10.85, (p = 3.55, A = 8.08, B = 0.00007, 
and C = —0.00037. The r.m.s. of the fit errors is 0.0176. The parameter values of 
the fit with (4) are = lO-Aug-2000, a = 0.64, u = 10.89, (pi = 3.42, 02 = 0.51, 
A = 8.08, B = 0.00026, C = -0.00080, and D = -0.00026. The r.m.s. of the fit 
errors is 0.0153. The formula including the harmonics reduces the r.m.s. of fit errors 
by 13.2%. 
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Fig. 24. The stock market index of Mexico (fine line) and its fits with the simple 
log-periodic formula (1) (thick line) and with formula (4) incorporating an harmon- 
ics (dashed line). The starting date for the fits is tstart = 04-Sep-2000. The param- 
eter values of the fit with (1) are = 08-Aug-2000, a = 0.80, uj = 10.11, c/) = 4.42, 
A = 8.70, B = 0.00032, and C = 0.00086. The r.m.s. of the fit errors is 0.0402. 
The parameter values of the fit with (4) are tc = 27-Jun-2000, a = 0.73, uj = 5.97, 
01 = 4.23, 02 = 1-65, A = 8.69, B = 0.00072, C = 0.00039, and D = -0.00115. 
The r.m.s. of the fit errors is 0.0391. The formula including the harmonics reduces 
the r.m.s. of fit errors by 2.6%. 
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Fig. 25. The stock market index of Indonesia (fine line) and its fits with the simple 
log-periodic formula (1) (thick line) and with formula (4) incorporating an har- 
monics (dashed line). The starting date for the fits is tgtart = 16-Sep-2000. The 
parameter values of the fit with (1) are tc = 09-Aug-2000, a = 1.06, to = 10.30, 
(j) = 5.83, A = 5.99, B = 0.00009, and C = -0.00021. The r.m.s. of the fit errors 
is 0.039. The parameter values of the fit with (4) are tc = 16-Nov-2000, a = 0.96, 
uj = 3.47, (pi = 1.33, (/)2 = 3.10, A = 6.04, B = -0.00025, C = -0.00050, and 
D = —0.00062. The r.m.s. of the fit errors is 0.0271. The formula including the 
harmonics reduces the r.m.s. of fit errors by 30.4%. 
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Fig. 26. The stock market index of Korea (fine line) and its fits with the simple 
log-periodic formula (1) (thick line) and with formula (4) incorporating an har- 
monics (dashed line). The starting date for the fits is tstart = 12-Oct-2000. The 
parameter values of the fit with (1) are tc = 15-Nov-2000, a = 1.37, to = 6.69, 
(j) = 5.67, A = 6.29, B = 0.00005, and C = -0.00005. The r.m.s. of the fit errors 
is 0.05. The parameter values of the fit with (4) are tc = lO-Aug-2000, a = 1.00, 
uj = 5.04, 01 = 1.90, (t>2 = 1.90, A = 6.24, B = 0.00046, C = -0.00023, and 
D = —0.00025. The r.m.s. of the fit errors is 0.0428. The formula including the 
harmonics reduces the r.m.s. of fit errors by 14.4%. 
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Fig. 27. The stock market index of Thailand (fine Une) and its fits with the simple 
log-periodic formula (1) (thick line) and with formula (4) incorporating an harmon- 
ics (dashed line). The starting date for the fits is tstart = 20-Sep-2000. The param- 
eter values of the fit with (1) are tc = 16-Aug-2000, a = 0.93, uj = 9.45, cp = 5.16, 
A = 5.61, B = 0.00058, and C = -0.00039. The r.m.s. of the fit errors is 0.0437. 
The parameter values of the fit with (4) are tc = 27-Jul-2000, a = 0.43, uj = 5.18, 
01 = 3.62, (j)2 = 2.22, A = 5.48, B = 0.02083, C = 0.00390, and D = 0.00662. The 
r.m.s. of the fit errors is 0.0354. The formula including the harmonics reduces the 
r.m.s. of fit errors by 19.0%. 
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Fig. 28. The stock market index of Russia (fine line) and its fits with the simple 
log-periodic formula (f ) (thick line) and with formula (4) incorporating an harmon- 
ics (dashed line). The starting date for the fits is tstart = Ol-Dec-2000. The param- 
eter values of the fit with (1) are tc = 08-Oct-2000, a = 0.92, u = 7.93, = 3.33, 
A = 7.24, B = 0.00276, and C = -0.00068. The r.m.s. of the fit errors is 0.0485. 
The parameter values of the fit with (4) are tc = 12-Oct-2000, a = 0.91, uj = 7.85, 
01 = 3.90, 02 = 0.17, A = 7.23, B = 0.00307, C = -0.00075, and D = 0.00010. 
The r.m.s. of the fit errors is 0.0457. The formula including the harmonics reduces 
the r.m.s. of fit errors by 5.7%. 
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Fig. 29. The French stock index, a bearish anti-bubble, fitted from tstart to Septem- 
ber 2002 for different choices of tstarti spanning from Jun-01-2000 to Dec-01-2000. 
One see that the fits are robust with respect to different starting date. 
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Fig. 30. The Australian stock index, a bearish anti-bubble, fitted from tstart 
to September 2002 for different choices of tstartj spanning from Jun-01-2000 to 
Dec-01-2000. One see that the fits are robust with respect to different starting 
date. 
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Fig. 31. The stock market indices of seven countries (US S&P 500, the Netherlands, 
France, Germany, Norway, UK, Spain) are superimposed by plotting the normalized 
values [p{t) — {p)]/cTp of each index as a function of time, where {p) is the mean 
whose substraction accounts for a country-specific translation in price and Up is the 
standard deviation for each index which accounts for a country-specific adjustment 
of scale. 
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Fig. 32. Cross-correlation coefficients between the USA S&;P500 index and four Eu- 
ropean stock market indexes (Netherlands, France, Germany and United Kingdom 
from top to bottom) for weekly returns on filtered prices in a moving three-month 
window. Surrogate tests obtained by reshuffling the returns show that the large 
peaks and troughs are statistically significant. One can observe a slow overall in- 
crease of the correlation coefficient C{t) over this decade. The variations are esti- 
mated to be 0.12 by bootstrapping simulations. 
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Fig. 33. Cross-correlation coefHcients between the USA S&;P500 index and four 

non-European indexes (Hong Kong, Australia, Russia and China from top to bot- 
tom) for weekly returns on filtered prices in a moving three-month window. While 
one can observe a slow overall increase of the correlation coefficient C{t) over this 
decade for Hong Kong and for Australia, no such trend is seen for Russia and 
China which are characterized by weaker correlations and more pronounced regime 
switches. The variations are estimated to be 0.12 by bootstrapping simulations. 
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Fig. 34. Monthly returns of four European market indices as a function of the 
monthly returns of the USA S&PSOO index. The black dots correspond to pe- 
riod [Jun-04-1997, Aug-09-2000] . The + correspond to the period [Aug-10-2000, 
Sep-04-2002]. 
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Fig. 35. Weekly returns of four non-European market indices as a function of 
the weekly returns of the USA S&PSOO index. The black dots correspond to pe- 
riod [Jun-04-1997, Aug-09-2000] . The + correspond to the period [Aug-10-2000, 
Sep-04-2002]. 
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Fig. 36. The synchronization factor R{t) between the USA S&PSOO index and four 
European indices (HL, FR, DE and UK), defined as the fraction of weeks with the 
same return signs over a moving 13-week window. The horizontal hues show the 
average of R{t). The variations are estimated to be 0.06 by bootstrapping simula- 
tions. 



64 



N»4 




1982 1984 1986 1988 1990 1992 1994 1996 1998 2000 2002 2004 



0.8 
0.6 
0.4^ 



-AU . . h J, At-A^iJ' 


,1 






1 r 


1 1 1 1 1 1 1 - 



1984 1986 1988 1990 1992 1994 1996 1998 2000 2002 2004 




2003 



2003 



Fig. 37. The synchronization factor R{t) between the USA S&PSOO index and four 
non-European indexes (HK, AU, RU and CN), defined as the fraction of weeks 
with the same return signs over a moving 13-week window. The horizontal hues 
show the average of R{t). The variations are estimated to be 0.06 by bootstrapping 
simulations. 
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Fig. 38. Synchronization index Q{t) (see text and expression (13) for its definition) 
between tlie USA S&P500 index and four European stock indexes (HL, FR, DE and 
UK), as a function of time in a running window of 65 trading days. The horizontal 
Unes show the average of Q{t) over the shown time interval. 



66 




67 



6000 



/= 0.841 < > CO = 5.28 




Log frequency 



